Faraday instability in a two-component Bose Einstein condensate 
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Motivated by recent experiments on Faraday waves in Bose Einstein condensates (BEC) we in- 
vestigate the dynamics of two component cigar shaped BEC subject to periodic modulation of the 
strength of the transverse confinement. It is shown that two coupled Mathieu equations govern the 
dynamics of the system. We found that the two component BEC in a phase mixed state is relatively 
more unstable towards pattern formation than the phase segregated state. 
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P_h ' Faraday waves are generated when the free surface of a fluid layer is subjected to a periodic vertical acceleration PJ. 
When the acceleration exceeds a threshold value, surface waves appear oscillating at half the forcing frequency. Re- 

■ cently, Faraday waves, which comes under the category of spontaneous pattern formation, a very general phenomenon 
studied in different fields of nonlinear science, were observed in a cigar shaped BEC by periodically modulating the 

r ^ 1, radial trap frequency Q. The radial modulation leads to a periodic modulation of the density of the cloud in time, 
which in turn leads to a periodic change in the non-linear interactions. This leads to the parametric excitation of lon- 
gitudinal soundlike waves (Faraday waves) in the direction of weak confinement. It has been shown theoretically that 
for a BEC, Faraday instability can be generated either by modulating the scattering length by Feshbach resonance 
Q or by modulating the trap frequency in the tight confinement direction (4)]. In both the cases, the dynamics are 
& ' governed by a Mathieu equation that is typical for parametrically driven systems. Floquet analysis reveals that a se- 
ries of resonances exists, consisting of a main resonance at half the driving frequency and higher resonance tongues at 
interger multiples of half the driving frequency . A theoretical analysis of Faraday instability based on a Mathieu-type 
analysis of the non-polynomial Schrodinger equation was given recently 0. From the perspective of phonon number 
occupation, the Faraday type modulation has been analyzed theoretically [6| . In general, the nonlinear spatiotemporal 
dynamics of BECs is attracting increasing interest in recent years, with a major focus on controlling structures like 
solitons by temporal modulation of the atomic scattering length [3]. The evolution of a BEC in a time dependent 
trap has been addressed earlier by some authors Q. 

The aim of the present paper is to analyze Faraday instability in a two-component cigar shaped BEC. The two 
components could be for example 87 Kb atoms in two different hyperfine states in different external trapping potentials. 
A periodic modulation of the radial trap frequencies would give rise to Faraday instability in both the components. 
The dynamics of this system is now governed by two coupled Mathieu equations. If £!i and f?2 are the natural 
frequencies of the two components then the parametric resonances will occur not only at 2f2i/n and 20,2 /n (n is an 
integer) but also at the combination frequencies (| ± Qi ± Q,2\/n). The behaviour of this system would then be similar 

■ to spontaneous pattern formation in miscible and immiscible binary classical liquids [9( . 

An important question is the stability of such Faraday waves in binary mixtures because depending on the atom- 
atom interaction strengths ga (intra species interactions) and , * j£i (inter species interaction) , binary mixtures can 
be phase separated (3111722 < 312) or phase mixed (311322 > 312) The nonlinear interaction can be conveniently 
manipulated by Feshbach resonances [ill ]. The purpose of the present paper is to analyze the influence of the 
interactions on the stability of the parametric excitations that are generated as a result of periodic modulation of the 
radial confinements of the two cigar shaped elongated components of the BEC held in a magnetic trap. In contrast 
to binary liquids, the interactions in the present system are tunable and one can in the same experiment crossover 
from one state to the other. In particular, we show that a phase mixed state is more unstable towards parametric 
excitations than a phase separated state. 
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II. THE MODEL 



We consider a situation in which we have a two component BEC (bosonic atoms of the same isotope but having 
different internal spin states, eg. 87 Kb atoms in hyperfine states F = 2,mp = 2 and F = l,mj? = —1). The starting 
point of our analysis is the coupled Gross Pitaeviskii (GP) equations for a two component confined BEC: 
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Figure 1: Resonance tongues of the parametric instability for the dissipationless phase separated state (gi — 0.5, 32 = 0.75). 
Shaded domains indicate where the two-component system is unstable. Starting from the lowest tongue, the resonances are 
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Figure 2: Resonance tongues of the parametric instability for the dissipationless phase mixed state (gi = 1.25, 32 = 1.5). 
Shaded domains indicate where the two-component system is unstable. Starting from the lowest tongue, the resonances are 

fil + ^2 

SI2, , Hi and 2^2, fii + SZa, 2fh. Clearly compared to Figure 1, nearby resonance tongues are overlapping and thereby 

the stable regions around the resonance tongues are reduced indicating that a lower value of the perturbation is sufficient to 
induce a Faraday instability in the phase mixed state as compared to a phase segregated state. 



lh ^tr = I + Vi(f,t) + g^N^l 2 + g 12 N 2 \^ 2 \ 2 - ^ Ui (1) 



dt y 2m 

lH ^f = {-^ 2 + V 2 (r,t)+g 22 N 2 \iP 2 \ 2 +g 21 N 1 \^ 1 \ 2 -112^2 (2) 
Here tpi, with i = 1,2 is the effective wavefunction of the i th condensate, with mass m and chemical potential 
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Hi. The interaction between the i condensate atoms is specified by gn = and that between 1 and 2 by 

2 m 

.912 = 921 = — ■• Here an is the intraspecies s-wave scattering length and ayi — a^\ is the interspecies scattering 

rn 

length. Ni is the number of particles of the i th species. Here aa(aij) > for a repulsive BEC, which we consider. 
Vi(r, t) is the trapping potential of the i th component. The wavefunctions are normalized according to / \ipi\ 2 dr = 1 
in the absence of damping. Here we are ignoring damping and by doing so, we do not loose the essential physics. 
Later qualitatively, we will discuss the role of damping on the stability chart and include it while studying the 
spatio-temporal dynamics of the BEC. 

When the condition .911522 < 3i2 is satisfied, the condensates are phase segregated. In the opposite limit i.e 
911921 > 9i2, the condensates are in a mixed state. Our study covers the ID case of a cigar shaped BEC extended 

777 

along the z direction, Vi(f, t) = — [Sl^ (t) (a; 2 + y 2 ) + f2 2 u .z 2 ], where Cl Wi << Vl ti is assumed which means that fl Wi 

and flti are the frequencies of the trap for the i th component along the weak and tight confinement directions, 
respectively. We assume that S7« i is subjected to periodic modulation: f2t 4 (t) = Qt 4 [1 + on cos fi it] , on << 1. Here 
f2j is the trap modulation frequency of the i th component. The radial modulation leads to a periodic change of the 
density of the cloud in time, which is equivalent to a change in the nonlinear interactions and speed of sound. In this 
paper we will always work either in the deep phase mixed regime or in the deep phase segregated regime. At the 
boundary separating the phase mixed regime and phase separated regime, the dynamics could be very complex as 
the trap modulations can continuously switch the system from one regime to the other. In terms of scaled variables, 



t = £!i2i, R = (X,Y,Z) = (x,y, z)/a±2, ^12 = 5- and ttj = a\i\] AnNiayi^i-, Eqns.(l) and (2) takes the following 
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dimensionless form 



lit = H v * + \< {x2 + y2) + \< z2 + I |mi ' 2 + |u2 ' 2 - "1 ( " 1 :::s s 



dt \ 2 " 2 12 v 7 2 W2 2 

n 2(^22 ^"Wi — /-, \ _ Mi 

where gi = , uj ti = —— , ui Wl = — — , Lu ti = w ti (l + on cosc^t), m = , W = 7^— • 

gi2 II12 "12 hitl2 "12 

As in [J|, we reduce the BEC 3D dynamics to an effective ID description using a multiple scale analysis pjj. We 
write the dimensionless wavefunction as: 

Uj - ^lexp {~^f{X 2 +F 2 )) fc^r),* = 1,2 (5) 
For a flat trap in a weakly confined space uj Wi = , consequently Eqns. (3) and (4) become 

2l ~g^ = \ ~~Q£1 + 5 , l' Zl *i 1^1 1 2 + <^7?, 1 02 | 2 - A > 01 +5l^t I «l COSCJir|0i| 2 0i (6) 



2t -g^ = \-QZ2 + 52^*2 |02| 2 + Wr|0i| 2 - ^ | '■ '2 ' .92^/ , <» 2 COS u-'o rj , ) 2 | <. - ( i \ 

i 

where w# = — — — and /i' = 2 f/Zj — uj t ■ ) , i — 1,2 

Wti + w t2 

In order to proceed analytically, we study Eqns. (6) and (7) for the simple case ol\ = ui = a, w\ = L02 — w and 
uj tl — u t . 2 — uj t . In this case lvr = (D t (l + acoswr). For a flat potential and in the absence of any perturbation 
(a = 0), (j>i = 0° is a constant independent of Z. Eqns. (6) and (7) yield: 

^H 1 -5iff2j 
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Figure 3: Farday pattern formation in space and time for the first condensate in the phase mixed state. The parameters are 
uj — 0.47 x 1CP 3 , a = O.f, 7 = 0.1 ,gi = 1.25 and 52 = 1.5. The time and lenght are in the units of Q12 and 012 respectively. 




Figure 4: Faraday pattern formation in space and time for the second condensate in the phase mixed state. The parameters 
are u> = 0.47 x 10 -3 , a — 0.1, 7 = 0.1 ,gi = 1.25 and g^ — 1.5. The time and lenght are in the units of Q12 and 0,12 respectively. 
Note that for the second condensate since the driving frequency is near one of its resonance tongues, the Faraday patterns 
generated are more compared to that in the first condensate. 



1-O12 Pi- 9102 x 

ml = — n T = <52 ( 9 ) 

wt{i--gig2) 

where /3j = 2(/2j — ui t ), i = 1,2. In the presence of the modulations, Eqns. (6) and (7) admits the following 
homogeneous states: 
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b-y = v Siexp ( — — sinwrj (10) 



/n 



'52exp ( — — sinwr ) (11) 



"I- 1 ~ff 2 I U0I2 _ I 1 -9i 



This yields, ^ = ,4 = 1 , |0° | 2 = | ^ | and | & , 

Now we wish to know whether the spatially homogeneous external driving fields are able to induce a spontaneous 
spatial symmetry breaking of this homogeneous state of the coupled BEC's. For that, we perform a linear stability 
analysis of Eqns. (6) and (7) by adding a small perturbation to Eqn. (10) in the form: 

4>i = $ (l + WiCosk Si Z) ,i = 1,2 (12) 

where Wi and k Si , (i = 1,2) is the complex valued amplitude and the wave-vector of the perturbation of the i th 
component respectively. Substituting Eqn. (11) into Eqns. (6) and (7) and linearizing with respect to Wi leads to the 
following coupled Mathieu equations for m — ReWi, (z = 1,2) 



d ui 
d 2 u 2 



(Q 1 + aQn coswr) u\ + Q12 (1 + acoswr) U2 = (13) 



(Sl| + CKQ22 cos wr) u 2 + Q21 (1 + acosLur)ui = (14) 



where Oi(k Si ) — k 2 t + 2<fe, i s t ne dispersion relation of the perturbation of the i th condensate in the absence 

of driving. Also q n = |^— — ^-|, q 22 = \y^- — ^r\, qu = qn/gi, <Z2i = 922/52, Q11 = k% q n /2, Q22 = k 2 g 2 2/2, 
{1-9x92) {1-9x92) 

Qx2 = ~k 2 Sl qx2/2,Q 2 x = ~k 2 S2 q2x/2. 

As known from the general theory of equations of the type (12) and (13), instability may occur for u> near twice 

the natural frequencies and their subharmonics, 2f2j/n., i = 1,2, n is an integer, and also close to the so-called 

combination frequencies and their subharmonics, | ± Cli ± Q,2\/n. The boundary curves for the instability regions 

must therefore emerge from these frequencies. The wavenumbers fc Si ,n corresponding to the resonance tongues near 



the natural frequencies and their subharmonics is fc Si . n = y — qu + \J qu + 4(nw/2) 2 . In order to find the stability 

diagrams for the two coupled Mathieu equations, we employ the method indicated in [13]. The method involves 
infinite determinants which are truncated and then the eigenvalues are found using MATHEMATICA, together with 
a condition on the characteristic exponent for the solutions on the boundary curves. The result for the phase separated 
case (51 = 0.5, 172 = 0.75) and the phase mixed case [g\ — 1.25, 32 = 1-5) is shown in Fig.l and Fig. 2 respectively. 
The shaded regions are the regions of instability. Clearly one notices that the two sets of resonance tongues, namely 

0,1,0,2, and 20i, 2^2,^1 + ^2 for the phase separated case are more widely spaced with larger stability 

region in between the tongues than the phase mixed case. The stable regions between the resonance tongues for the 
phase separated mixture, extends as far as upto a — 0.2 while those for the phase mixed case extends only upto 
a — 0.1. If the detuning of the driving frequency is of the same order for two adjacent modes then these modes can 
enter in competition. The competition between nearly degenerate modes has been shown to lead to chaotic behaviour 
in fluids Since the resonance tongues are more closely spaced in the phase mixed state, the mode competition 
will be stronger for this case. 

UJ2_U0|2 

Now, we study the time and space evolution of the deviation from the initial density fi = — — » % — by including 

\9i\ 

phenomenologically a damping term proportional to 7 — - in Eqns. (13) and (14). Here 7 is the damping coefficient. 

The coupled Mathieu equations (13) and (14) are solved using MATHEMATICA. Figures (3) and (4) illustrates the 
time and space evolution of f\ (density deviation of the first condensate) and f% (density deviation of the second 
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condensate) in the phase mixed state respectively. The parameters are u> = 0.47 x 10 -3 , a = 0.1, 7 = 0.1 , g\ = 1.25 
and <?2 = 1-5. The periodicity of Faraday patterns generated is always more for the condensate with lower natural 
frequency. We found that this observation is also true for the phase separated case. The influence of interactions 



on the periodicity of pattern formation is clear from the expression fc Si . n = y —qu + \J qu + 4(mj /2) 2 . Tuning the 



interactions, controls the k Si and hence the periodicity. Also noticed is the fact that onset of pattern formation is 
not immediate but takes a certain time due to damping. The fact that the phase mixed state is relatively more 
unstable compared to the phase segregated state was also revealed by the spatio-temporal dynamics (figures not 
shown). We found that for a given driving lo and perturbation a, the growth rate of the perturbations in the phase 
mixed state was much higher than that in the phase separated state. For long enough times in the presence of a 
continuous modulation, the excitations will grow and eventually the condensates will be destroyed. This contrasts 
with the everlasting periodic revivals of the spatial modulations in one component ID case [1 51 ] . The spontaneous 
pattern formation for combinations frequencies is much more complex and simple analytical formulations cannot be 
done and would require a full numerical simulation of the coupled equations (1) and (2). 



In conclusion, we have demonstrated that the dynamics of Faraday instability in a two component elongated cigar 
shaped BEC is described by two coupled Mathieu equations. The dynamics depend on the two body intra and inter 
species interactions. In contrast to the single component case, we found that parametric resonances can occur not 
only at integral multiples of half the natural frequencies of the two components but also at integral multiples of half 
the combination frequencies. The stability chart in the {uj, a} plane reveals that the a phase separated binary mixture 
is more stable towards pattern formation than the phase mixed mixture. This was confirmed by the spatio-temporal 
dynamics of the BEC where we found that the initial perturbation grows in a phase mixed state more rapidly than 
in the phase segregated case. 
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